Superconducting Spiral Phase in the two-dimensional t-J model 
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( ^^ ) We analyse the t — t' — t" — J model, relevant to the superconducting cuprates. By using chiral 

perturbation theory we have determined the ground state to be a spiral for small doping 8 <C 1 
near half filling. In this limit the solution does not contain any uncontrolled approximations. We 
evaluate the spin-wave Green's functions and address the issue of stability of the spiral state, leading 
to the phase diagram of the model. At t' = t" = the spiral state is unstable towards a local 
enhancement of the spiral pitch, and the nature of the true ground state remains unclear. However, 
, for values of t' and t" corresponding to real cuprates the (1,0) spiral state is stabilized by quantum 

£SJ ■ fluctuations ("order from disorder" effect). We show that at 8 ~ 0.119 the spiral is commensurate 

with the lattice with a period of 8 lattice spacings. It is also demonstrated that spin-wave mediated 
superconductivity develops in the spiral state and a lower limit for the superconducting gap is 
. derived. Even though one cannot classify the gap symmetry according to the lattice representations 

(s,p,d,...) since the symmetry of the lattice is spontaneously broken by the spiral, the gap always 
has lines of nodes along the (1, ±1) directions. 



C3 



O 



T3 

o 



X 



PACS numbers: 11.30.Er, 31.30.Jv, 31.30.Gs 



I. INTRODUCTION 



C 

The t — J model has been suggested to describe the essential low-energy physics of the high-T c cupratesi 2 ^. 
Formation of spirals in the doped t — J model was proposed by Shraiman and Siggiai They showed that the pitch of 
the spiral is proportional to the hole concentration 6, see also Ref£. The idea that the spiral state is the ground state 
Cn| . of the doped Heisenberg antiferromagnet attracted the attention of theorists£i£*2'i2*ii*i2*i^, however the question 
of stability of the state remained controversial. According to Refs. 6 ' 8 - 10 the spiral state is unstable toward a local 
enhancement of the spiral pitch. On the other hand the analysis of Ref. 9 indicated that the spiral state is stable. 
A semiclassical analysis of stability of the (1,0) spiral state was performed in Ref£. According to this analysis the 
state is marginal (zero stiffness), which effectively indicates that the state is unstable. A complete stability analysis of 
spiral states in the Hubbard model was performed in Ref— also using the semiclassical approximation. According to 
their analysis in the leading in powers of doping approximation the (1,1) spiral is always unstable and the (1,0) spiral 
is always marginal in agreement with*. Recently the interest in the spiral state was renewed because of the strong 
experimental indications that at small doping the cuprates behave as spin glasses or even exhibit some kind of magnetic 
ordering. A summary of the available experimental data on one of the superconducting cuprates, La2_5Sr^Cu04, is 
given in Ref— . The data also show that magnetic ordering and superconductivity coexist. The spin glass behavior is 
consistent with the spiral scenario: since doping is not uniform the pitch of the spiral is varying from point to point 
and hence on large scales it leads to spin glass behavior—. 

In the present work we analyse the stability of spiral states within the RPA approximation. The approximation is 
parametrically justified for S <C 1. In this part of the work on the technical side we follow the approach developed 
by Igarashi and Fuldc 9 . However, contrary to them and in agreement with RefsA&, we conclude that in the "pure" 
t — J model (i! = t" = 0) the spiral state is unstable with respect to a local enhancement of the spiral pitch. We find 
that relatively small values of t' and t" stabilize the spiral order, with a uniform hole distribution. Our results are 



then consistent with numerical (DMRG) results indicating that inhomogeneous (striped) phases disappear with the 



increase of further-neighbor hoppingsiSiii. In our approach the instability of the spiral is closely related to the fact 
that the hole dispersion is almost degenerate along the face of the magnetic Brillouin zone. As soon as the degeneracy 
is sufficiently lifted by t' and t" the instability disappears, leading to a stable (1,0) spiral state. Within the semiclasical 
approximation the effective stiffness of this state is zero in agreement with£*Ai. However spin quantum fluctuations give 
rise to a nonzero positive stiffness and hence favor stability via an "order from disorder" mechanism. For parameter 
values corresponding to real cuprates, the pitch of the spiral is proportional to doping, and at S ~ 0.119 the spiral 
becomes commensurate with the lattice with period 8, in agreement with the experimental data of Tranquada et 

am. 

The possibility of s- and d-wave pairing between mobile holes due to a static distortion of the Neel background 
was pointed out in Ref—. Such a distortion leads to an infinite set of very shallow two-hole bound statesSi. In 
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RefiS^ it was shown, under the assumption that Neel order is preserved under doping, that in the t — J model there is 
spin- wave mediated superconducting pairing in all waves except the s-wave. The pairing was found to be maximum in 
the d-wave channel. A numerical calculation performed in Ref.2 3 . under the same assumption showed that the pairing 
was quite strong. However, the picture of pairing 22-23 was not consistent because the staring point of the analysis 
was the (unstable) Neel background. In the present paper, after proving the stability of the spiral state, we consider 
superconducting pairing on this state. Using an approach similar to that of Ref^ we find a superconducting pairing 
instability and show that the gap always has lines of nodes along (1, ±1) directions. We estimate analytically the lower 
limit for the superconducting gap. One cannot classify the gap according to the lattice symmetry representations 
since the symmetry is spontaneously broken by the spiral. 

In essence the method we use is chiral perturbation theory which allows the treatment of strong interactions and 
is exact in the long wave-length limit 24. The small parameter of the approach is doping near half filling, S <C 1. We 
cannot determine reliably whether S — 0.15, S = 0.1 or 6 = 0.05 is sufficiently small to justify our calculations, but we 
claim that at sufficiently small S only the long range dynamics is important and the approach is parametrically justified. 
The analysis of pairing within the chiral perturbation theory by construction complies with the Adler relation^ and 
hence with the argument presented by Schrieffe:»2& concerning pairing via exchange of Goldstone excitations. 

The structure of the paper is as follows. In Section [H] we calculate the single hole dispersion and the quasiparticle 
residue for different values of t' and t" , To do so we use the self consistent Born approximation. After that we 
perform the RPA analysis and demonstrate the instability of the Neel order, at arbitrary small hole concentrations, 
towards formation of spirals. In Section III we consider the hole dispersion in different spiral states and compare the 
total energies of the states. The RPA analysis of stability of the (1,0) spiral state is performed in Section IV, and in 
Section V we calculate the reduction of the spiral on-site magnetization due to doping. Section VI is devoted to the 
superconducting pairing in spiral states. We present our conclusions in Section VII. 



II. SINGLE HOLE PROPERTIES AND RPA PROOF OF INSTABILITY OF THE NEEL ORDER UPON 

DOPING 

The Hamiltonian of the t — t' — t" — J model is: 

(ij)a (ijl)<? (ih)cr (ij)cr 

c\ a is the creation operator of an electron with spin a (a ="[", j) at site i of the two-dimensional square lattice, (ij) 
represents nearest neighbor sites, (iji) - next nearest neighbors (diagonal), and {1J2) represents next next nearest sites. 
The spin operator is Si = \c\ a a a pCip, and the number density operator is rii = ^ CT c^ ia Ci G . The size of the exchange 
measured in two magnon Raman scatterin g 27 ' 28 ! is J = 125meV. Calculations of the hopping matrix elements have 
been performed by Andersen et a&. They consider a two-plane situation and the effective matrix elements are 
slightly different for symmetric and antisymmetric combinations of orbitals between planes. After averaging over 
these combinations we obtain: t — 386meU, i! — — 105meU, t" = 86meV. From now on we set J = 1. In these units 
we have: 

t = 3.1, t' = -0.8, t" = 0.7 (2) 

An analysis of angle-resolved-photoemission spectra for the insulating copper oxide S^CuC^C^ performed in Refi^ - 
with the Hamiltonian CQ),© shows an excellent agreement with experiment for both the single-hole dispersion and 
for the photoemission intensity. This analysis is based on the Self-consistent Born Approximation (SCBA) 31 - 32 . The 
approximation works very well due to the absence of single loop corrections to the hole-spin- wave vertesSi24. In the 
present calculation of single-hole properties we follow the approach of Refi^i and use the SCBA. 

It is well known that without doping (half filling) the Hamiltonian is equivalent to the 2D Heisenberg model 
which has antiferromagnetic (Neel) order. There are two sublattices: sublattice a with spin up and sublattice b with 
spin down. Hence a hole created in the system has a pseudospin index a or b. The bare hole operator di is defined so 
that d\ oc Cij on the a sublattice and cx Cii on the b sublattice. In the momentum representation: 

d ^ = ^ Nil/ l +m) Y. c ^> ^ = ^(l/2 2 + m) g C ^' (3) 

where N is the number of sites and m = |(0|Sj Z |0)| « 0.3 is the average sublattice magnetization. The quasi- 
momentum k is limited to be inside the magnetic Brillouin zone: 7k = ^(cos k x + cos k y ) > 0. The index a (b) can be 
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considered as pseudospin. Rotational invariance is spontaneously broken, but nevertheless the pseudospin gives the 
correct value of the spin z-projection: a corresponds to S z = — 1/2 and b corresponds to S z = +1/2. The coefficients 
in J2J provide the correct normalization: 

(0K.4J0) = JV(1/2 2 +m) E<Olc! T c it 10) = v ±- { 0\l + S tM \0) = 1. (4) 

For spin excitations the usual linear spin- wave theory is used (see, e.g. the review paper—). Spin- wave excitations 
are described by operators a q and /3_ creating spin waves with S z — —1 and S z — +1 respectively. The momentum 
q is restricted inside the magnetic Brillouin zone. The operators are defined by the equations: 

|]TS+e-^ « u^ + v^, (5) 

Jl^Sre^ « Vq < + Uq /5_ q . 

The spin- wave dispersion and the parameters of the Bogoliubov transformation diagonalizing the spin- wave Hamilto- 
nian are: 



'£4 



1 1 



Wq 2' 



vq = -signed) \ 

Hopping to a nearest neighbor site in the Hamiltonian (JIJ leads to an interaction of the hole with the spin-waves: 

H h>sw = ^2 9k,q (4+qa d k;,aq + 4+qfAa/?q + H.C^j , (7) 



k.q 

with the vertex g k ,q given b y 33 i 34 



5k,q = (0|a q rf k 6|fft|4+qal°) =4 *Y^(TkWq + 7k+qWq). ( 8 ) 

The vertex (?k,q is independent of t', t" because these parameters correspond to hopping inside one sublattice. In 
agreement with Adler relation^ the vertex vanishes at q — + 0. We calculate the hole Green's function using the SCBA 
approximation^ii 3 ^ 3 ^ 3 ^. This gives the quasiparticle dispersion ek and the quasiparticlc residue Z^. In the vicinity 
of the dispersion minima, ko = (±7r/2, ±7r/2), the quasiparticle residue is: 

Z k «Z = Z ko , (9) 

and the dispersion can be approximated as 

2 2 2 2 

e k w const + fa-fi + /3 2 ( 7k f » const + + 2 ^ ^ /?iy + ft y ■ (10) 

Here 7^ = ^(cosfc x cosfcy), p — k ko, the component pi is orthogonal to the face of the Brillouin zone and the 
component p 2 is parallel to the face, see FigQ] The results of our calculations at t = 3.1 and for the second neighbor 
hoppings within the intervals — 1 < t 1 < and < t" < 1 can be fitted as: 

t = 3.1, 

f3i = 1.96 + 1.1U' + 0.06i' 2 + 2.7CM" + 0.53i" 2 + 0.50i'i" , 

#2 = 0.30 - 1.33t' - 0.19i' 2 + 2.80t" + 1.06i" 2 - O.Ut't" , 

Z = 0.29 + 0.055i' + 0.195i" . (11) 
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These formulae agree with the results of Refsi2&2i at t' = t" = 0, as well as with Refi^S for t', t" given by J5J. 
The quasiparticle-spin-wave interaction differs from the bare interaction Q only by the presence of the quasiparticle 
residues and reads: 

H qp . sw = ^ (hl+^hkbOq + frk+qfc /l ka/3q + H.C?j , (12) 
k,q 

where /i ka and are quasihole creation operators, and the vertex M q has the form (we set N = 1 for convenience): 

M q = v/Z k Z k+q5k , q « -2 7 ^Zt^- . (13) 

Here qi is the component of q orthogonal to the face of the Brillouin zone, and we assume that k ss ko and q -C 1. 
It is convenient to introduce the fields 7r q and A q instead of a q and /3 q : 

7T q Q! q /3_ q , 

A q = a q + /3l q . (14) 
The Green's functions of these fields are defined as: 



C„Kq) - -i / die ia "(T[7r q (t)7r q (0)]) , 

/•+OC 

D Xn (u,q) = -i die iwt (T[A q (t)A+(0)]) . (15) 



In the absence of interactions they are: 

I>W ( Wl q) = D% (u,, q) = „ a J^ + . Q ■ (16) 

The index n in l|15|l . i|l()|l labels them as normal Green's functions. According to Eq. l|12|l the field 7r q interacts with 
the quasiholes: 

Hgp,n = A/ q (^k+qa^kb7r q + ff.C.) , (17) 
k,q 

while the field A q remains idle with respect to this interaction. The interaction (|17|) generates a loop correction to the 
spin-wave Green's function shown in Fig. [21 The holes are fermions and occupy four half pockets or two full pockets 
as shown in Fig[3] At a given momentum there are two states with different pseudospin, so the Fermi energy ep and 
the total Fermi motion energy Ep (per lattice site) read: 



E F = j^/JmS 2 , (18) 
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FIG. 2: TYq-spin-wave polarization operator 
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FIG. 3: Hole pockets and Fermi surface for the Neel state and for the spiral state with Q oc (1, 0). 



where 8 -C 1 is the hole concentration. The calculation of the polarization operator Fig. [21 is straightforward and 
gives: 



p{, , _ ac n \ - \^ m 2 \ ^ 2 ^ e P ~ £ p+q) n p( 1 ~ n p+q) _ 2 5 / 2 Z 2 f 2 1 ^ 2 pff > 



pockets 



p+q 



(19) 



pockets 



Pfcq) = ^ [t 2 -2i? ev /(t 2 /2 + ^) 2 -2£Ft 2 ] , F(0,0) = 1 



Here i 2 = /3ig 2 + /?292 an d n p i s the Fermi distribution function. Re stays for real part. Summation over two full 
pockets shown on the right side of Fig. [3] has to be performed. Note that what is qi for one pocket is q2 for the other, 
and vice versa. At lo = and q -C pf ~ V<5 the polarization operator (|19fl is particularly simple: 



, . 2 5 / 2 Z 2 t 2 
P(0,q) = j== q 

71" V P1P2 



(20) 



The polarization operator is always proportional to q, as a direct consequence of Goldstone's theorem. A crucial 
point is that the polarization operator (|20|l is independent of i5. Taking into account the one- loop diagram Fig. [21 the 
7r-magnon Green's function (|15f) becomes: 



A™ (w, q) 



2uj„ 



o\ - 2w q P(w, q) + iQ 



(21) 



The stability of the system requires that all poles of the Green's function l|21|l are at positive uj 2 . Using Eqs. i|2U|) 
and i|ll|) one can see that this criterion is violated for q <C Pf ~ V<5 at arbitrary small doping, because (cf. with 
Refsi 1 - 36 ) 



cj q + 2w q P(w = 0,q) = 2q' 1- 



2 j- 2 



< 



(22) 



This result signals an instability towards formation of spirals. 
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In the presented RPA proof of the Neel state instability we have used the following approximations: (1.) The 
single hole properties (dispersion and quasiparticle residue) have been calculated without account of other holes. (2.) 
The incoherent part of the hole Green's function has been neglected. (3.) Interactions between quasiholes have been 
neglected. These approximations are parametrically justified at sufficiently small 6 since (|2U|1 is independent of S. We 
observe that doping influences the behavior of the system only at momenta q ~ \/6 -C 1. At the same time all integrals 
in the SCBA are convergent at large momenta, q ~ 1. Consequently doping gives a negligible correction (oc S) to the 
dispersion. The same is true for the incoherent part of the hole Green's function and for the hole-hole interaction. 
The leading correction from these effects is oc 6 In 5, see Refi^i, and hence it is also negligible for sufficiently small S. 



III. SPIRAL STATES, HOLE DISPERSION AND TOTAL ENERGY 

In a spiral state there are still two sublattices, sublattice a and sublattice 6, but the spin at every site j of each 
sublattice is rotated by an angle Bf 

\tpj) = e^ n - CT/2 | t) if je a sublattice , 
\iPj) = e ie i n -^ 2 \l) if j£ b sublattice, 
6j = Q • rj . (23) 

Here Q <C 1 is the vector of the spiral, a is the vector of Pauli matrices, and n is an axis of rotation. The direction of 
n can be arbitrary in the plane orthogonal to the spin | f) in Eq. I|23l) . n = [n\, Ti2, 0) = (cos a, sin a, 0). To first order 
in Q the small rotation l|23|) does not influence the spin- wave dispersion and the hole-spin- wave interaction considered 
in the previous section. The only effect which appears at first order in Q is the possibility for a hole to hop between 
nearest neighbor sites of the lattice. Using Q and i|23|) one can easily find the Hamiltonian describing this hopping: 

H Q = -t Y, (^) [<e- fa 4 f4i6 d ia - ie ia d\ a d i+ ^ b ] . (24) 

Here £ is a lattice vector directed from a given site i S a to the nearest neighbor i + £ G b. One can rewrite (|24H in 
momentum representation using the quasihole operators h^: 

Hq = -t y~](Q x sin k x + Q y sin k y ) (e~ la d^d^ + e ia d\, a d*A 

k 

-> -Z k tY(Q*sink x + Q y smky)(e- ia hl b h ka + e ia hl a hu b ) . (25) 

k 

The quasiparticle residue is exactly the same as the one found in the previous section by the SCBA. Indeed, the 
effects considered in the present section are relevant to the long range dynamics, so the corresponding momenta are 
very small, Q oc 5 <C 1. On the other hand all the integrals in the SCBA are convergent at large momenta and 
therefore not sensitive to the long range dynamics. The axes x and y coincide with the crystal axes. Note that the 
axes 1 and 2 used to define the phase a are not related to the crystal axes. The Hamiltonian (|25|l mixes the states 
hl. a and h^. Therefore the effective Hamiltonian matrix for a given k takes the form: 

tt _ ( £k -Z k te~ la [Q x smk x +Q y smk y ]\ , . 

eff ~ \-Z k te la [Q x smk x + Q y smk y } e k J ' [Z0> 

where ek is found in the previous section. In the vicinity of the points ko = (±7r/2, ±tt/2) one can use Eqs. @ and 
HI U| l and hence H e f / takes the following form: 

H <"-\y/2Zte<«Q x e k ) ' (27) 

where Q\ is the component of Q orthogonal to the corresponding face of the magnetic Brillouin zone. We denote by 
■0 k and the creation operators that diagonalize this Hamiltonian 

1$ — — (h^ - p^^h 
^k — ^ V ko kb 

1 

V2 



= -^K a + e-^hl b ) , (28) 
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where e tfi = y^e 1 "- The corresponding dispersions are: 



e k - =ek- V2Zt|Qi| , 
£k+ = e k + \/2Zi|Q 1 | . 



(29) 



We will see that \Qi\ is always large enough, so only the states corresponding to the e k _ branch (^-band) are filled by 
holes. According to Eq. (|29|l the total kinetic energy gain due to the spiral is AEkin = —V2Zt\Qi\S per lattice site. 
On the other hand the spiral increases the magnetic energy. The variation of the magnetic energy per lattice site is 
AE magn = ^p s Q 2 , where p s = Z p /A w 0.18 is the spin stiffness, and Z p w 0.72 is the renormalization factor due to 
higher l/S'-corrections 38 . Finally, the Fermi motion of holes Ep oc S 2 , also contribute to the total energy. Altogether 
the total energy of the spiral state (per site) with respect to the undoped antiferromagnet is: 



E= ^ Ps Q 2 - V^ZtlQilS + E F 



(30) 



There are two possibilities to minimize the total energy: for a spiral directed along a diagonal of the lattice Q oc (1, 1), 
and for a spiral directed along a crystal axis of the lattice Q oc (1,0). For the diagonal spiral the effective Fermi 
surface consists of two half pockets or one whole pocket, see Fig. The hole states do not have a pseudospin index 
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FIG. 4: Hole pockets and Fermi surface for the spiral state with Q oc (1, 1). 
and hence the Fermi energy ep and the total Fermi motion energy Ep (per site) in this case (Q oc (1, 1)) read: 



E F = Tt^/frfoS 2 ■ 



(31) 



For the spiral along a crystal axis the effective Fermi surface consists of four half pockets or two whole pockets, 
similarly to the Neel case, see Fig. [21 However the hole states in a spiral do not carry a pseudospin index and hence 
the Fermi energy and the total Fermi motion energy Ep in this case (Q oc (1,0)) are: 



f-p = TTy/J^JhS 



E F = -VM^ 2 



(32) 



Using Eqs. (|30|l . (|31() . and l|32l) we find the following expressions for the energy of the spiral state relative to the 
energy of the undoped antiferromagnet: 



Qoc(l,l): E {hl) = X - Ps Q 2 - V2ZtQ5 + 7r^fafoS 2 
Qoc(l,0): E {m = Ip^-ZtQS+^-y^S 2 . 



(33) 
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Minimizing E with respect to Q one finds: 



Q 0,(1,1): Q = ^ls, £ (1 , 1) = (-^! + 7 r v /^) ( 5 2 

Ps I Ps J 

Q oc (1,0) : Q = —6, E (m = I — - — + -Vftft U 2 . 

Ps I 2p s 2 J 



(34) 



Eqs. (fTTT^f) and (|3"4*)l agree with Ref£, see also Refi^. We see that -E(i.i) = 2^(1,0), so if these energies are positive 
then the state (1,0) has lower energy than the (1,1) state. It should be noted that the energy -£(1,0) is always 
lower than the energy of the doped Neel state (i.e. the difference between i^i.o) and Ea. H18(l is always negative). 
When £'(i l i) is negative then this state has lower energy. However, in this case the corresponding compressibility 

— Qfri'^ is also negative which implies an instability towards a state with an inhomogeneous charge distribution. 
Consequently the (1,1) spiral state is always unstable in agreement with the conclusions of RefsASiii. On the other 
hand the comressibility of the (1,0) state is always positive (in agreement with Ref>ii). Using the results of our SCBA 
calculations, summarized in Eqs. (|11|1 . as well as Eqs. (|34|) . one can easily find the region where the state (1,0) is 
realised. Our results are summarized in Fig. The energy £(1,0) vanishes along the solid line, and the spiral state 
is unstable in the right bottom corner of the phase diagram. To complete the analysis of stability of the (1,0) spiral 
state we must calculate the spin- wave polarization operator ("transverse stiffness"). We leave this analysis to the next 
section. 




FIG. 5: Phase diagram of the t — t'—t" — J model at t = 3.1 and small doping, S <C 1. The bottom right corner corresponds to 
the pure t — J model, t' — t" = 0. The top left corner corresponds to parameters from Refr—, Eq. The region of stability 

of the superconducting (1,0) spiral phase is shown. 

For values of t, t' and t" from Refi^. 9 , Eq. @, the (1, 0) spiral is stable. At this point Z = 0.38, Eq.^TJ; an d hence 
according to Eq. the magnitude of the spiral vector is: 

Q = 6.585 , att = 3.1, t' = -0.8, t" = 0.7 (35) 

The spiral is commensurate with the lattice if Qn = ir, where n is an integer number. The corresponding hole 
concentration 5 n — 7r/(6.58n). This results in an effective antiferromagnetic structure with period Ar. The period 
is Ar = n for odd values of n, and Ar = 2n for even values of n. In Table 1 we present the values of S n and the 
corresponding periods of the commensurate antiferromagnetic structure for several values of n. 
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n 


3 


4 


5 


6 




0.159 


0.119 


0.095 


0.079 


Ar 


3 


8 


5 


12 



Table I. Values of the hole concentration S n at which a commensurate antiferromagnetic structure with period Ar is 

established. 



Notice that for 5 n = 0.119 the period is 8, in agreement with experimental data 18 . The spin structure at this doping 
value is shown in Fig. [5] In section VII we discuss further the possible connections of our results to the physics of 
the superconducting cuprates. 

t \ /_\ t 

\ t I 

FIG. 6: (1,0) Spiral state for 8 = 0.119, corresponding to Q = n/A. 



IV. SPIN- WAVE GREEN'S FUNCTION IN THE SPIRAL STATE AND RPA STABILITY ANALYSIS 

OF THE STATE. 



The stability analysis of the (1,0) spiral state requires a calculation of the spin-wave Green's functions l|15l) . similar 
to the calculation performed in Section[n]for the Neel state. The spin wave A q , Eq. I|14l) . remains unchanged since it 
does not interact with the holes. The interaction of the spin wave 7r q with the holes is given by Eq. I|17|) . but one has 
to take into account the fact that the correct quasiparticle states are represented by the operators tpk and ip-^ instead 
of hk a and /ikb, Eq. \2H\i. The diagrams shown in Fig. contribute to the 7r-spin-wave polarization operators. The 



p+q 




p+q 




p+q 




-q 



p+q 




FIG. 7: Contributions to ix -spin-wave polarization operators. A and C describe transitions within the ip-band, and B and D 
describe transitions between the ip- and tp-bands. 



diagrams Fig. 0A and B contribute to the normal polarization operator. The filled state is always ifa, however the 
excited state in the loop can be either ipk+q or </?k+q, see Eq. I)28[l. providing the difference between the diagrams Fig. 
13 A and B. Since pseudospin is not conserved there is also an anomalous polarization operator. The contributions to 
the anomalous operator are given by the diagrams Fig. {7\C and D. The anomalous polarization operator gives rise to 
the anomalous 7r-spin-wave Green's function defined as: 

/+oo 
die iwt (T[7rl q (t)7r q (0)]) . (36) 
-OO 
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A straightforward calculation of the polarization operators Pa and Pc, taking into account transitions within the 
same subband, gives (we set uj = i£ from now on): 

r>i c \ 2iu D ct- \ 2(e p — e p+q )?ip(l — rt p+q ) \/2Z 2 t 2 1 ^— ^ i-vic \ toi\ 

PA^) = -e -Pc(Cq) ^ E TL-^ + ^—^j L ^,q). (37) 



pockets 



pockets 



The function P(£,q) is defined in Eq. (|T§|| . Eq. i|3T|) differs from (fT^f) by the coefficient 1/4 only. The summation 
has to be performed over two full pockets of the (1, 0) spiral state, Fig. We remind the reader that q± is the 
component orthogonal to the face of the Brillouin zone. Next we evaluate the polarization operators Pb and Pjj 
describing transitions between subbands: 



Mi 



P B «,q)=e 2 ^/W,q) = E -f E 



2(e 



p e p+q 



A)n F 



V2Z 2 t 2 1 



* ^.JL^^- (38) 



q) = i [* 2 + 2A - 2i? ev /(t 2 /2 + A + i0 2 -2e F t : 



*(0,0) = §1 



where A = 2\/2Zt\Qi \ is the splitting between the subbands, see Eq. I|29l) . Taking into account Eq. (|34() for the (1,0) 
state we obtain: 



72 f 2 

According to Fig. [7|and Eqs. (|37|) . I|38l) the normal and the anomalous polarization operators are: 

V2ZH 2 1 



(39) 



P«W,q) - P A {iti,<l) + P B {it,<l) 

Pa(^,q) = fc(^,q) + -PD(^,q) 



pockets 

V2Z 2 t 2 e~ 2 ^ 



kVpip2 q 



E ^H^q) + q)] 



(40) 



pockets 



The corresponding Dyson's equations for the Green's functions read: 

Amfa) - D$ (?) + D$ (q)P n (q)D nn (q) + (q)P a (q)D na (q) 
D na (q) = D^{-q)P n {-q)D va {q) + D^{-q)P*{q)D m {q) , 



(41) 



with Dml given by Eq. . By solving these equations we obtain the normal and the anomalous Green's functions 
of the 7r-field: 



Am(w = i£,q) = 



-2a; q [C 2 +^ + 2cj q P Tt ^,q)] 
[C 2 + w q + 2a, q P n (i£,q)] 2 - 4c q |P Q ^,q)r 

2u; q P a *(^,q) 

[e + Wq + 2^ q P„ (it, q)] 2 - 4 Wq | P Q q) | 5 



(42) 



As already discussed in section |HJ a necessary condition for the stability of the system is the absence of poles in the 
Green's functions at negative uj 2 (positive £ 2 ). The most dangerous is the regime of very low frequencies and momenta 
(w <C ep, q <C pf ~ VS). Therefore the criterion of stability is: 



[^ 2 +2 Wq P l (0,q)] 2 >4^ 2 |P a (0,q)| 2 
Using Eq. (|40|l the stability criterion can be rewritten as: 



1 > 



2ZH 



2,-2 



T^VJhJhps 



zn 2 



2Z 2 i 2 



Z 2 t 2 



(43) 



(44) 



According to Eq. I|34l) the expression 1 — tt\/ fiifoPs IZ 2 ^ 2 is negative in the (1,0) phase since P(i.o) is positive. 
Therefore the criterion (|44|l reads: 



1 - 4p s = 0.28 > , 



(45) 
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meaning that the spiral phase is always stable. It is important to note that the phase is stable only due to spin quantum 
fluctuations, and the stability is an order from disorder effect. Without account of the fluctuations (semiclassical 
approximation), p s = 1/4, and hence the effective transverse stiffness vanishes, 1 — 4p s = 0. In this case the system 
becomes marginal, or in essence unstable, as pointed out in Refs. 5,11 . Thus the spin quantum fluctuations make 
the (1,0) spiral phase stable. We stress that the spin quantum fluctuation effects in an undoped antifcrromagnet are 
known quite accurately from numerous approaches (spin- wave, series expansions, Monte-Carlo, etc.). Within the chiral 
perturbation theory, used in the present work, one can express the leading terms in the hole density via the parameters 
of the undoped system. In this sense we perform an "exact" account of the quantum fluctuations. However, we do not 
calculate terms subleading in the hole density because such a calculation cannot be performed without uncontrolled 
mean field approximations. It is possible that at sufficiently large doping the spiral phase becomes unstable due to 
such subleading terms. One possible type of instability is towards formation of a noncoplanar state, as suggested inAi. 
We cannot determine the exact value of doping S which is small enough to justify our calculations, but we claim that 
at sufficiently small 8 our approach is parametrically justified, and consequently Fig. [5] represents the phase diagram 
of the t — t' — t" — J model. We emphasize that in this regime the (1,0) spiral phase is stable for values of parameters 
corresponding to real cuprates. 



V. ON-SITE MAGNETIZATION IN THE SPIRAL STATE 



In the undoped t-J model the spins order in a staggered collinear pattern and the value of the staggered magneti- 
zation is: 

M = \ (S Z ) \ w 0.303 . (46) 

In a spiral state the on-site magnetization follows a spiral pattern and besides that the value of the magnetization is 
reduced compared to l|46|) . We calculate now the reduction of the magnetization. Using the spin- wave operators (J5J 
and neglecting corrections proportional to Q 2 oc 8 2 one can rewrite the on-site magnetization in the following form: 



M 



\{S z a -Si) = 0.303 - 2]T— ([c4a q +4/3 q -7q(«qi9-q + at^-q)]> 

= 1-2V- <[4a q + /9-di8i q -7q(«q/3-q + «i/9lq)]> • (47) 

q W 

For small doping the deviation of M from the value f4l)|l is due to quantum fluctuations at small momenta, g < 1. 
Hence one can replace 7 q — * 1 and rewrite (|47|l in terms of 7r-field averages and then the normal Green's function 
Eqs. CSJ, (X3 : 

1-2V— ( 7r t 7rq ) = l-2iV— D m (t = -0,q) . (48) 



M 



Finally, using the explicit expression (|42ll for the Green's function we find the reduction of the on-site magnetization 
in a spiral: 



M = 0.303 




-u£ + 2w q P n (i£,q) 



AM = -4 > / ^ { ' q V - —-^ . (49) 

' ' °" ' ' 1 2o, q P n (i£,q)] 2 -4c q |P a (^,q)| 2 £ 2 + ^ J 

A straightforward numerical integration of this equation with parameters corresponding to Eq. J2J gives the on-site 
magnetization plotted in Fig. |H1 

For hole concentration S = 0.12, corresponding to a commensurate spiral (Fig. |SJ), the value of the magne- 
tization is M « 0.06. It is quite close to the reported experimental values in Lai.88Sro.i2Cu04 4 ° as well as 
Lai.48Ndo.4Sro.i2Cu0^iS. The overall decrease of M as a function of doping is also consistent with experiments^. We 
observe that the spiral state disappears at S c ~ 0.16, beyond which point (<5 > S c ) spin rotational invariance is restored. 
The existence of a magnetic quantum critical point inside the superconducting region is an issue of considerable cur- 
rent interest (the spiral states support superconductivity as discussed in Section VI). It has also been suggested 
that the phase emerging upon destruction of non-collinear magnetic order could be of RVB type and exhibit electron 
fractionalizationi2i2£, although in the present work we can not address the structure of the magnetically disordered 
region S > 6 C . 
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FIG. 8: On-site magnetization in the (1,0) spiral versus doping for t — 3.1, t' — —0.8, and t" = 0.7. 

Using the expression for the polarization operator (|40|l one can show that AM can be expanded in powers of 5 in 
the following way: 

AM = -aSln(l/6) + bS + cS 2 + ... . (50) 

This equation shows that the expansion is nonanalytic in S. We have calculated the coefficients a and b (as well as c) 
numerically for several values of t' and t" , and tabulated them in Table II. 

Only the logarithmic term in the expansion l|5U|l is parametrically justified within the chiral perturbation theory 
which we use in the present work. This is due to the fact that the integral in Eq. (|49|l behaves as 6 J dq/q and 
consequently the "b"-term in l|50(l is related to the upper limit of the integration, i.e. to the short distance dynamics 
(<J~tt). 



spiral state 


t' 


t" 


a 


b 


(1,0) 


-0.8 


0.7 


1.3 


-1.9 


(1,0) 


-0.45 


0.45 


1.2 


-1.3 


(1,0) 


-0.3 


0.25 


1.0 


-0.7 



Table II. Coefficients a and b in the S-expansion \5(]\) for the reduction of the on-site magnetization for different 

values of the hopping parameters t' and t" . 

Eq. 149(1 overestimates the contribution of large momenta because it does not take into account the reduction of the 
quasiparticle residue away from ko = (±7r/2, ±7r/2). On the other hand it underestimates this contribution since it 
does not take into account the incoherent part of the hole Green's function. To estimate the uncertainty related to the 
short distance dynamics we have calculated the reduction of the nearest- neighbor sites spin-spin correlator (Si • Sj) 
in the same way, i.e. assuming no variation of the quasiparticle residue and without the incoherent contribution. The 
result of this calculation agrees within 30% with the available numerical data^. This indicates that the values of the 
coefficient "b" in Table II are somewhat reliable. We stress once again that the leading "a" -term in 1|50|) is due to the 
long distance dynamics only and therefore it is parametrically justified and reliable. 

VI. SUPERCONDUCTING PAIRING IN THE SPIRAL DUE TO SPIN- WAVE EXCHANGE 

Now we investigate the possibility of superconducting pairing due to the spin-wave exchange shown in Fig. We 
consider pairing between the states ijjp, Eq- H28fl . inside one full pocket (see the right hand side in Fig. [Jlfor the (1,0) 
state). The corresponding many-body wave function has the form: 

i*> = n ( u p + v p^- P ) i°) - ( 5i ) 
p 

and represents pairing of spinless fermions. Within the full pocket description the typical momentum transfer in the 
diagrams Fig. [§] is q = p — p', q = |q| ~ pp oc \f5 <C 1. In the description with half pockets (left hand side in 
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Xq qy 
s x » >'' ^ 

-p -p' -p -p' 

FIG. 9: Spin-wave exchange between holes leading to superconducting pairing 



Fig. 13) the momentum transfer is close to the antiferromagnetic vector q w G = (±tt, ±tt). We use the full pocket 
description which is more convenient, and in fact Eq. 1)51(1 already assumes such a description with V_ p = — Vp. 
Note that this condition does not mean that we consider negative parity pairing. Since parity is defined in the full 
magnetic Brillouin zone, in order to consider parity one needs to transfer back from the full-pocket description to 
the half-pocket description. Such a transition is related to translation by the antiferromagnetic vector G. The wave 



function changes sign under such translation, = — Vk + Q. Therefore V- 



-V p implies that V_ k = +V-y and the 



pairing has positive parity. For a discussion of this peculiar symmetry property we refer the reader to RefsiSi* 2 ^. 

Thus in the pairing channel the typical momentum transfer is q ~ pp oc y5. At such momentum transfer the 
spin- wave Green's function ((42(1 is substantially different from the bare one l)16[l. Its spectral weight ImD w (u)) is 
shifted from u = aj q = \/2q oc y/S down to lj ~ 6f oc 5. Unfortunately in this regime we cannot find an analytical 
solution for the pairing. Even at small S one needs to solve the Eliashberg equations numerically and this is outside of 
the scope of the present work. However we can put analytically a lower limit on the pairing and study its symmetry. 
To do so we replace the renormalized Green's function (|42|l in the diagrams Fig. El by the bare Green's function JT^J. 
We call this approximation the "bare approximation" - it is similar to the one used in Rcf. 2 — . The typical energy 
transfer in the diagrams of Fig. [5] is u ~ t F oc 5 <C w q oc \/6. Therefore when calculating the diagrams Fig. one 
can set u> = in the spin- wave Green's function 1(16(1 . This means that in the "bare approximation" it is sufficient to 
consider only the static interaction. By calculating the diagrams Fig. [5] with account of Eqs (|16f) . (|17fl . and l|28() we 
find the effective pairing potential due to the spin-wave exchange: 



(dir) 



M 2 

q 



72 .2 <?1 



= -8ZH*Z± = -8Z 2 t 2< p- 
T (P- 



P'i) 2 



P') 2 



(52) 



Recall that pi is the component orthogonal to the face of the magnetic Brillouin zone. Note that in agreement with 
the general Adler relation 2 ^ M q — > 0, q — > and hence the pairing potential remains finite at q = 0, as was pointed out 
by Schrieffer— . In addition to the diagrams in Fig. there are also exchange diagrams which differ by permutation 
of the legs. The total effective pairing potential due to the spin-wave exchange is then: 



Vp, p / = V {dir) - v {exch) = 8Z 2 t 2 



(Pi-P'i) 2 | (Pi+P'i) 2 
(p-p') 2 (p + p') s 



-A_ p reads: 



The BCS equation for the superconducting gap A p 



p' 



p' 



(dir) Ap' 



(53) 



(54) 



where E p — (e p — fi) 2 + A 2 . Note that using the interaction 1)53(1 one has to put the coefficient 1/4 instead of 1/2 

in the BCS equation to avoid double counting in the wave function (|51|l . Alternatively one can use the standard form 
with 1/2 and with the direct interaction ((52(1 only. In the weak-coupling limit the gap is small compared to the Fermi 
energy, therefore Eq. ((54(1 can be solved in the vicinity of the Fermi surface (Pxpi/2 + fapli/'Z ~ £f) with logarithmic 
accuracy. The solution is discussed in Refi 2 ^ and we directly present the result: 

A p = A sc sinm^ , 



A SC = Ce F e~ 1/g ™ 
AZ 2 t 2 

9rn 



yJMfh 



Here A$c is the maximum value of the superconducting gap. The angle cf> is defined as: 



(55) 
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FIG. 10: d-wave like symmetry of the superconducting gap in the (1,0) spiral state. 



to = 1, 2, 3, .. is an integer number, and C ~ 1 is a constant. Eq. (|55J) represents a family of solutions, however the 
pairing is maximum in the channel with m = 1. In this case, near the Fermi surface, the gap has the simple form 
(with pi and p2 defined in Fig. [Q: 

y/(fil/p2)Pi+P2 

We observe that there is a line of nodes along the (1, 1) direction. 

The solution (|55J) only describes the pairing within one full pocket while for a (1,0) spiral (which is of primary 
interest in the present work) there are two full pockets (Fig. I2J). In this case one needs to solve the BCS equation 
numerically (i.e. take into account the pocket-pocket scattering) in order to determine the symmetry of the solution. 
This has been done for a similar problem in Refi^ and the expected result is shown in Fig. 1 101 where we also have gone 
back to the full Brillouin zone (half pocket) description. Since the symmetry of the lattice is spontaneously broken 
by the spiral, the "d-wave" classification has lost meaning. In principle this should be reflected as an asymmetry of 
the pockets themselves; this asymmetry is small for 6 <C 1 and is not shown in Fig. (|10fl . 

The pairing l|55|) obtained in the "bare approximation" is rather weak. As we have already mentioned the use of 
the exact spin-wave Green's function Ij42(l is expected to enhance the pairing substantially due to the spectral weight 
shift towards low frequencies; we plan to discuss the full numerical solution of the Eliashberg equations in a future 
work. However there exists an additional contribution to pairing (not accounted for in the solution iJSSJ) which can 
be estimated already in the weak-coupling limit. This involves pairing via the upper ip-band (ghost band) and can be 
included by replacing the many-body wave function (|51|l by: 

p p' 

To account for this effect we calculate the second order correction to the spin- wave mediated interaction llol'li between 
the ■(/'-fermions. The correction is given by the diagram in Fig. 1111 The dot in this diagram is given by Fig. [5] where 




FIG. 11: Second order correction to the pairing interaction between the tp-fermions due to virtual excitations to the empty 
ip-band. The dot represents the spin-wave exchange shown in Fig. 

the "flavor" of the fermion is changed from ip to ip in each spin- wave vertex. Calculating the diagrams in Fig. |5|with 
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account of Eqs. (|16fl . (|17fl . and (|28|l we find that the dot is described by the same Eq. (|52|l but with a sign minus: 
"dot" — Af 2 /u> q = 8Z 2 £ 2 (p! — li) 2 /(p — l) 2 . Hence the second order correction Fig. ITT1 reads: 

Wp, P ' - -( 8Z * ) (p _ 1)2 2ei + 2A-2 ei , (p'-l) 2 ' (59) 

where ei is given by Eq. jlU|>. and and A are given by Eqs. I|32p. I|39l) for the (1,0) spiral state. The evaluation of 
the sum in Eq. (|59|) is particularly simple for the isotropic case, /3\ = fa. We will consider this case for the purpose 
of an estimate of the effect. Using the expansion with Chebyshev polynomials T n (x): 

i-^^( 1 + 2 g^) ' ™ 
and keeping only the first angular harmonic which is important for pairing in the m = 1 channel, we obtain: 

(Pl - Zl) 2 1 V< I , . , N / C1 S 

-«»(^ + ^). ( 61 ) 

where p< = min(p, I) and p> = max(p, I). Using this representation and having in mind that in Eq. <|59|) \p\ = |p'| = 
Pi? we perform the integration in l|59|) and find the correction to the pairing interaction: 

x f dir) (8Z 2 t 2 ) 2 ( A-c F , A tp A \ . . 

<SU o' = — ,o a 1 ln a !~ a m — x ( cos ^p cos <?V + sm <^p sm <?V ) • ( 62 ) 

v,v 62/Kp V e_F A — £f A — £f/ 

We remind the reader that the above formula is valid in the isotropic case, = j3\ = 02- The correction (|62(l gives 
rise to an approximately 10% enhancement of the superconducting coupling constant g\ , Eq. (|55|l : 

ff * ^^A^ A + e^AX (63) 



51 27T/3 \ A — £f A — 

This enhancement is an interesting qualitative effect related to the "ghost" band, but numerically its influence on 
pairing is relatively weak. Taking into account both Eq. Q55JI and Eq. we obtain the following estimate for the 
pairing in the dominant, m = 1 channel (we use /3i = 2.9, /?2 = 3.8 corresponding to the hoppings from Eq. J2J), and 
also reinstate the Heisenberg exchange J): 

A sc ~ W- 2 SJ. (64) 

For ten percent doping this weak-coupling formula produces a T c cx Asc of the order of several Kelvin which is 
about 10 times smaller than transition temperatures in real compounds. The most important pairing enhancement is 
expected from the renormalization of the spin-wave Green's function, to be discussed in a separate work. 



VII. CONCLUSIONS AND DISCUSSION 



We have studied the phase diagram of the t — t' — t" — J model close to half-filling. To determine the single hole 
properties (one hole injected in an antifcrromagnetic background) we use the self consistent Born approximation. 
This approximation is not parametrically justified for the t — J model, but it is known that it works remarkably 
well for the single hole properties. The crucial observation is that these properties are not related to the long range 
dynamics: all integrals in momentum space are convergent at large momenta (quantum fluctuations at distances 1-2 
lattice spacings). Finite doping brings nontrivial long-range dynamics into the problem. To determine the dynamics 
we use the chiral perturbation theory which involves an expansion in powers of doping near half filling, 8 <C 1. The 
efficiency of the chiral perturbation theory is in essence a consequence of the dimensionality of the problem (2+1 
D) and Goldstone's theorem. We certainly cannot determine the exact value of doping so that it is small enough to 
justify the approach, but we claim that at sufficiently small 5 the approach is parametrically justified. We show that 
the Neel state is unstable with respect to decay to spiral states as soon as doping is introduced. The analysis of the 
stability of the spiral phases is performed within the RPA approximation, and all non-RPA corrections are suppressed 
by powers of 5. 

We find that at t' = t" = the spiral state is unstable toward a local enhancement of the spiral pitch, and 
consequently at that point the nature of the true ground state remains unclear. However we have shown that for 
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values of t' and t" corresponding to real cuprates the (1,0) spiral state is stable. The phase diagram of the model at 
t = 3.1 (we set J = 1) is shown in Fig. |SJ For hole concentration S w 0.119 the (1,0) spiral is commensurate with 
the lattice with a period of 8 lattice spacings, in agreement with experimental data. Fig. |H1 shows the dependence of 
the spiral on-site magnetization on doping. We demonstrate that spin-wave mediated superconductivity is developed 
above the spiral state and derive analytically a lower limit for the superconducting gap. One cannot classify the gap 
according to the lattice representations Ai, A2, Bi, B2, and E ("s","d","p",...) since the symmetry of the lattice is 
spontaneously broken by the spiral. However, the gap always has lines of nodes and a symmetry similar to d-wave. 

Finally we briefly discuss the possible connection of our results to experiments in the cuprates. Magnetic order 
has been observed in the superconducting states of both L^-a-^NdaSrgCuO ^i 19 and L^-^SraCuC^ 14 . This order 
seems to be particularly enhanced near the hole concentration 6 = 1/8 where a commensurate structure with a period 
of 8 lattice spacings is observed. In addition, static charge order (with a period 4) has been observed at least in one 
of the materials 1 ^*— This has lead to proposals that the state near 6 = 1/8 is a collinear spin density wave (SDW) 
coexisting with a charge density wave (CDW) at twice the wave vector—. Such an interpretation of experiments rules 
out a non collinear spiral configuration of the type shown in Fig. El since the density of holes is expected to be uniform 
in the spiral case. On the other hand the question whether (static) charge order is a generic feature of the cuprates 
is far from being resolved and experiments in many cuprates are interpreted as showing fluctuating (i.e. dynamic) 
order in the charge sectoo 4 ^ For example the existence of (static) charge order in La2-5Sr«5Cu04 is far less cleani^ 
although it also exhibits commensurate magnetic peaks around 5 = 1/8. 

In the light of the above we find our result that the spiral (1,0) state is commensurate with the lattice (period 8, 
Fig. ©) for 6 = 0.119 very promising. It seems to be consistent with the data in the magnetic sector (although 
the interpretation of Ref— rules out a spiral state in that material). Moreover, it is known that doping introduces 
disorder, in turn leading to glassy features in the magnetism; at least at low doping (in the normal state) those 
features can be explained well within the spiral scenario— In our opinion it would be also very interesting to study 
the density response of the system in the spiral state. In the context of the t — J model the density response has 
been recently investigated in the Neel stated, where fluctuations at energies of order J were found. In a spiral state, 
due to the presence of the ghost (empty) band, we expect that even lower energy charge density fluctuations will 
exist, although a perfectly static CDW order seems impossible to stabilize. It is certainly more likely that some 
kind of order in the charge sector appears (if anywhere at all) around the point t' = t" = 0, where the spiral state 
is inherently unstable. If we indeed interpret the region marked as "unstable" on our phase diagram Fig. [S] as a 
candidate for such order, then the effect of increasing t' and/or t" is to drive the system towards the (stable) spiral 
order, with a homogeneous charge distribution. This occurs for \t'\/t > 0.18 (at t" — 0). From this point of view 
our results are similar, at least superficially, to the DMRG results describing the destruction of stiped phases by 
second-neighbor hoppingi§*ii, where a critical value of t' necessary to destabilize the stripes was found, and is quite 
close to our estimate above. 
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